Dynamics and nonequilibrium states in the 
Hamiltonian mean-field model: A closer look 
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We critically revisit the evidence for the existence of quasistationary states in the globally coupled 
XY (or Hamiltonian mean-field) model. A slow-relaxation regime at long times is clearly revealed by 
numerical realizations of the model, but no traces of quasistationarity are found during the earlier 
stages of the evolution. We point out the nonergodic properties of this system in the short-time 
range, which makes a standard statistical description unsuitable. New aspects of the evolution 
during the nonergodic regime, and of the energy distribution function in the final approach to 
equilibrium, are disclosed. 
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The globally coupled version of the XY model, given 
by the Hamiltonian 



■H = K + V 
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[l-cos(0.-0,)], (1) 



has recently attracted considerable attention ||, |[ ^, ||, 
^, 0, ^, ^ . This model, usually referred to as the Hamil- 
tonian mean-field (HMF) model, describes an ensemble 
of N rotators characterized by their angles 9i and con- 
jugate momenta pi. In contrast to the usual XY model, 
interactions affect all pairs of rotators in the ensemble. 
The HMF model has been studied in connection with the 
emergence of collective self-organized states dynam- 
ical behavior near phase transitions Q|, and anomalous 
diffusion in phase space among other dynamical and 
thermodynamical aspects. 

The canonical thermodynamical description of the 
HMF model can be completely carried out [|l], ^ . In this 
description, it is useful to introduce the "magnetization" 
M = (cos 0„ sine,) = {M,,My) = Mexp(i$) 

and the average energy per particle U = E/N, with 
E = Ti.. The canonical approach predicts that the values 
of the magnetization and the average energy at equilib- 



rium, Mcq and Ucq, are related by 



T 1 



(2) 



where T is the canonical temperature. The system under- 
goes a ferromagnetic-like, second-order phase transition 
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at Tc = 1/2, from a state with M^q 7^ at low tem- 
peratures to a state of vanishing magnetization at high 
temperatures. The energy at the transition is C/c = 3/4. 
These predictions are mostly confirmed by the numerical 
solution of the HMF equations of motion at fixed total 
energy, i.e. in a microcanonical scenario. In this situa- 
tion the equilibrium temperature is defined through the 
kinetic energy 



T — K 



(3) 
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and U = Toq/2 -f (1 - M^^)/2 [cf. Eq. §)]. 

Deviations between the canonical predictions and the 
microcanonical (numerical) results were however re- 
ported in the region of energies just below the critical 
point, 0.5 < U < Uc, where numerical realizations re- 
veal extremely slow relaxation toward the asymptotic 
state ||. Such deviations, observed when the en- 
semble of rotators is prepared in specially chosen ini- 
tial conditions — namely, the so called water-bag initial 
conditions — have been associated with the existence of 
long-lived states where, in the thermodynamical limit 
N —* 00, the system would spend asymptotically long 
times. In that limit, canonical equilibrium would never 
be reached. It has been suggested that the HMF long- 
lived states — which would replace canonical equilibrium 
in the thermodynamical limit — would be well described 
by the equilibrium distributions predicted by Tsallis's 
generalization of thermodynamics |^, This gen- 
eralization, in fact, relaxes the assumption of extensiv- 
ity of the Boltzmann-Gibbs formulation, and is therefore 
presumably expected to describe the equilibrium statis- 
tics of systems with long-range interactions , as in the 
Hamiltonian (|l|). 

Water-bag (WB) initial conditions fix 0^ = for all 
i, and the momenta pi are chosen at random from a 
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uniform distribution in such a way that J^iPi — ^ ^^'^ 
K — X^iP?/-^ ~ NU. For these initial conditions, in 
fact, the potential energy vanishes. For a typical numer- 
ical realization at energies just below the critical point 
Uc, the dynamics of the HMF model can be summarized 
as follows. In the first stage, whose duration is essen- 
tially independent of the system size N and equals a few 
time units 0|, there is a rapid broadening of the dis- 
tribution in 9. This implies an abrupt growth of the 
potential energy from its initial value V{t = 0) = 0, 
to a value close to that predicted at canonical equilib- 
rium, Vcq = N{1 — Mgq)/2. The kinetic energy drops 
accordingly, from K{t = 0) = NU to a value close to 
Kcq — NToq/2, and the initial unbalance between K and 
V is drastically reduced. As explained below, the details 
of the evolution in this stage of kinetic relaxation depend 
strongly on the specific realization of the WB initial con- 
dition. Typically, however, the potential energy exceeds 
Vcq and attains a maximum, while K reaches a mini- 
mum. From then on, the system enters a stage of slow 
relaxation toward canonical equilibrium. 
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FIG. 1: Nonequilibrium temperature T{t) = 2K(t)/N of the 
HMF model for U = 0.69 and N = 500 (diamonds), 10^ 
(squares), and 10'' (circles), from Fig. 1(b) of Ref. ^ Note that 
data are plotted against t — tg, with to = 100. The horizontal 
line indicates the canonical equilibrium temperature Tcq = 
0.476. The inset shows the same data against t. The scales 
of the axis in the inset are the same as in the main plot. 

The slow-relaxation regime of the HMF model was re- 
ported in the first systematic study of this system In 
Ref. ^ it was quantitatively characterized by analyzing 
the relaxation of the nonequilibrium entropy S{t) toward 
its equilibrium value Scq as a function of N. The differ- 
ence \S(t) — Scq\ was shown to approach an asymptotic 
value of order N~^^^ within a relaxation time ~ N. 
The same scaling has recently been reported for the re- 
laxation time of the nonequilibrium temperature, T{t) = 
2K{t)/N [0. The states visited by the system during the 



slow-relaxation regime have been named quasistation- 
ary states, presumably by analogy with the quasiequi- 
libria reported for one-dimensional gravitational systems 
|l2| , 13 and similar long-lived states in many-body sys- 
tems with long-range interactions Jl^ . The name sug- 
gests that the time scales associated with the evolution 
of such states are much longer than those of any other 
regime during the whole evolution. However, the evo- 
lution of S{t) toward S^q reveals no dynamical regime 
other than slow relaxation |l^ . 

Latora, Rapisarda and Tsallis (LRT), on the other 
hand, have claimed to demonstrate the entity of quasita- 
tionary states (QSS) from the analysis of the evolution 
of the nonequilibrium temperature T{t) Q. For clarity 
in the discussion, we reproduce in the main plot of Fig. |l| 
their figure 1(b) of Ref. ^ where T{t) is plotted against 
time for several values of A^, as obtained from numerical 
integration of the HMF equations of motion for U = 0.69 
|l6| . This value of the energy corresponds to the largest 
observed deviation of numerical results with respect to 
the canonical equilibrium prediction ||l|, In the plot, 
the initial stage of kinetic relaxation is disregarded. For 
short times, we observe a wide plateau where T{t) is ap- 
parently constant, but substantially differs from its equi- 
librium value. At long times, T{t) has reached Tcq and re- 
mains at equilibrium. These two plateaus are connected 
by a crossover stage where T{t) exhibits an inflexion point 
and its variation is seemingly faster. Consequently, LRT 
conclude that the slow-relaxation regime consists in turn 
of a first stage where the system is "trapped" in a QSS, 
and a subsequent crossover toward canonical equilibrium. 
Since the position tc of the crossover shifts to larger 
times as A^ grows, tc ~ A^, LRT deduce that in the ther- 
modynamical limit the QSS will be observed at all times, 
and that canonical equilibrium will never be attained. As 
we argue in the following, however, the presence of the 
plateaus and, thus, the existence of a QSS in these spe- 
cific realizations of the HMF model, is an artifact of the 
peculiar way chosen to plot the data in Fig. ^ 

In the first place, in order to exclude the initial stage 
of kinetic relaxation, LRT choose to plot T{t) as a func- 
tion oi t — to, with to = 100; instead of simply cutting 
off the time axis at to. Such procedure, which would be 
innocuous in a linear time scale, has far-reaching conse- 
quences in the logarithmic scale of Fig. |l]. In fact, the 
point t = to becomes shifted to — oo in the logarithmic 
variable, which considerably enhances the impression of 
having a plateau at small t. In the inset of Fig. |^ we 
plot T{t) against its original variable — also in logarith- 
mic scale — showing that the short-time plateau shrinks 
by one decade. The existence of a crossover state join- 
ing the two plateaus is itself an illusion created by the 
logarithmic scale, even in the original variable t. Indeed, 
it can be easily shown that the graph of a function with 
monotonic first derivative can display an (otherwise inex- 
istent) inflection point by the simple expedient of using 
a logarithmic scale in the horizontal axis jl^ . 

In order to avoid these undesirable effects and to re- 
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veal the true nature of the evolution of T{t) in the slow- 
relaxation regime, we opt for using linear time scales. 
Figure || shows, plotted with diamonds, the same data 
as in Fig. for N = 500 up to i = 2500, while the in- 
set shows the whole temporal range. We see that there 
are no traces of a short-time plateau, except perhaps for 
a small interval around t = 500. It might however be 
argued that the vertical scale in this plot has been exag- 
geratedly amplified. To decide over this point, we have 
used a Pade-like approximation to fit the data over the 
whole time range, shown as a curve in the inset. The 
quality of the fitting is quite acceptable. Its extension to 
the short-time range, also shown as a curve in the main 
plot, proves that in this range the data exhibit the same 
trend as for larger times. The evolution in the whole 
range is therefore uniform, and we find no arguments to 
attribute the condition of quasistationarity to the short- 
time interval. As for the small interval around t — 500, 
we also include in Fig. || a data set from Rcf. ^ calcu- 
lated with higher numerical precision (open circles) [ p^ . 
These data reveal a minimum in T(t) at i w 300 and 
confirm the absence of any plateau. Figure || shows the 
data from Fig. |]for iV = 10^ and 10^ in linear time scale. 
Again, we find no evidence of the existence of plateaus 
or QSS. 
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FIG. 2: The same data as in Fig. § for iV = 500 (diamonds), 
with linear time scale. The inset shows the whole time range. 
The curve is a Pade-like approximation over the entire range. 
Open circles are higher-precision results from Ref. ^ 

Let us mention that, more indirectly, the observation 
of anomalous diffusion in the slow-relaxation stage of the 
HMF model [|| confir ms that no special regimes exist 
within that stage. In fact, superdiffusion in the unfolded 
0-space, 9 G (— oo,-foo), which is characterized by an 
anomalous dispersion law (A6'^) ~ t", is recorded in 
the whole slow-relaxation stage with a constant exponent 
a ~ 1.4. The transition to ordinary diffusion (a = 1) oc- 
curs at the same time as the system reaches canonical 



equilibrium. The short-time plateau of Fig. [| and the 
crossover to canonical equilibrium are not detected by 
the transport properties of this system. 
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FIG. 3; The same data as in Fig. |i| for (a) N -- 
N — 10*, with linear time scale. The inset in 



10^ and (b) 
a.) shows the 

whole time range, while the curve is a Pade-like approxima- 
tion over the entire range. 



The data shown in Figs. ^ and || exhibit an inconsis- 
tency regarding the behavior of the nonequilibrium tem- 
perature in the first stage of the slow-relaxation regime. 
In fact, while in the results for TV = 500 (circles) and 
for N — 10^ we observe a well-defined minimum in T{t), 
low-precision data for N = 500 (diamonds) and results 
for N = 10^ do not display such minimum. To clar- 
ify this aspect of the evolution, we have performed ex- 
tensive numerical calculations for the same value of the 
energy, U = 0.69, maintaining the relative error in the 
total energy conservation around /S.E/E ^ 10"'^, which 
makes it possible to average over many realizations of 
the water-bag initial conditions (5 x 10^ for N — 500 
to 10'^ for N = 3000). Our results for the nonequilib- 
rium temperature T(t) are shown in Fig. 0, as a function 
of the rescaled time t/N . Each dot has been obtained 
as an average in time, over O.IA^ to 0.2A^ time units. 
The arrows indicate, for N = 1000 and 3000, the typical 
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dispersion in the values of T{t) after averaging over the 
whole set of realizations. Fluctuations are indeed large, 
as expected near a phase transition and due to the highly 
chaotic HMF dynamics ^, which makes it difficult to 
obtain reliable numerical results. In any case, our nu- 
merical calculations consistently confirm the existence of 
a minimum in the nonequilibrium temperature, and the 
subsequent steady growth of T{t). The inset shows the 
position of the minimum, t,nin, as a function of N, in- 
cluding also the value for N — 10^ taken from the data 
of Fig. ^(b). The straight line corresponds to a power-law 
fitting, tmin (X N", which yields i/ = 1.7 ± 0.1. Though 
the plot spans less than two decades in N, the fact that 
^miii grows faster than linearly with the system size is 
clear. 

Meanwhile, no plateau is observed, at least, of the kind 
shown in Fig. |l|. A plateau may of course be compul- 
sively ascribed to the zone around the minimum of T(t) 
at ^min: but it would be far from coinciding in nature 
with the QSS reported by LRT. In fact, although the 
minimum becomes broader for larger N, it also shifts 
to longer times, which hardly insures that a stationary 
(nonequilibrium) temperature would be observed in the 
limit N oo. 
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FIG. 4: Nonequilibrium temperature, T{t) = 2K{t)/N, for 
U = 0.69 and three values of A*', as a function of the rescaled 
time t/N. The arrows indicate the typical dispersion of T{t) 
after averaging over a large set of realizations (~ 10^ to ~ 10^, 
see text). For — 500 the dispersion is negligible in this 
scale. The inset shows the time of the minimum in T{t) as a 
function of N. The straight line is a power-law fit, with slope 
1.7. 

The well-defined minimum of the nonequilibrium tem- 
perature T{t) at tmin establishes a quite natural bound- 
ary between short-time and long-time dynamical ranges. 
It is interesting to analyze in more detail the short-time 
evolution, with emphasis in the angle and momentum 
distribution of the ensemble in /i-space, where at each 



time the system is represented by a set of N points in 
the {9,p) plane. One-particle distributions of both an- 
gles and momenta have received particular attention in 
previous work 0, ||, |^, 0, so that such representation is 
useful for comparison. Figure ^ shows four snapshots, 
at different times, of the /i-space distribution in a single 
(typical) realization of the HMF model with TV = 10"*^ 
rotators, at U = 0.69. For clarity, the origin of the an- 
gle variable, 9 — 0, has been chosen to coincide in each 
plot with the phase <& of the total magnetization M. In 
this way, the center of the distribution coincides approx- 
imately with the center of the plot. Take into account, 
however, that $ performs, for short times, diffusive-like 
motion Q. The WB initial condition would here corre- 
spond to a uniform distribution along a vertical segment 
at = 0, with \p\ < \/6U ~ 2.03. As advanced above, 
we observe a rapid broadening of the distribution in the 
angle coordinate 9, accompanied by a moderate collapse 
in p. At t — 100, much of the ensemble has become 
strongly mixed in /i-space, due to the highly chaotic na- 
ture of the dynamics for this value of U ]l|, ||]. There 
is however a substantial fraction of the ensemble, repre- 
sented by the zones where points remain ordered, where 
the strong correlations imposed by the WB initial con- 
dition persist. This fraction is mainly situated near the 
center of the distribution, in a region of relatively small 
momentum, but ordered arrays of points are also seen in 
other zones of /i-space. For these times, thus, the system 
can be thought of as a mixture of two phases, namely, a 
strongly mixed, "gaseous" phase and a highly correlated, 
"condensed" phase. 
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FIG. 5: Snapshots in /t-space at four different times of a sin- 
gle realization of the HMF model with water-bag initial con- 
ditions (A'' = 10*, U — 0.69). In the four plots the horizontal 
and vertical axes span the intervals {—n,n) and (—2,2), re- 
spectively. 

The distribution of the condensed phase in /i-space re- 
sults to strongly depend on the specific realization of the 



WB initial condition. This sensibility to the initial con- 
ditions, illustrated in Fig. |[ is again a consequence of 
the chaotic dynamics that, as time elapses, amplifies the 
originally small variations between different realizations. 
An important byproduct of this property, combined with 
the persistence of the initial correlations, is that the con- 
densed phase is not ergodic, in the sense that averages 
over realizations of the WB initial condition yield a poor 
description of the statistical properties over time of any 
single realization. The lack of ergodicity of the condensed 
phase makes in turn the whole system noncrgodic in this 
short-time range. An illustration of this overall lack of 
ergodicity is provided by the short-time evolution of the 
nonequilibrium temperature. In any single realization, 
after its initial drop, T{t) displays chaotic oscillations of 
appreciable amplitude. Due to the sensibility to initial 
conditions, however, the phases of such oscillations differ 
between realizations, and the oscillations disappear upon 
averaging. This observation should warn of the averaging 
procedures used in previous work to study statistical 
properties such as the nonequilibrium temperature and 
the momentum distribution function in the initial regime 
of kinetic relaxation. The warning applies also to the av- 
erage data of T{t) shown in our Fig. ^ at the shortest 
times. Before the nonergodic condensate disappears, in 
fact, a (standard) statistical approach to the HMF model 
results to be of limited applicability. 



times, we have defined the condensate as formed by the 
rotators i such that the distance to i -I- 1 in /x-space. 
Si = is lower than a given 

threshold S. We fix the threshold as the average distance 
between contiguous pairs at t = 0, S = 2\J 6U /TV. 
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FIG. 6: Snapshots in /i-space at time f = 30 of four real- 
izations of the HMF model with water-bag initial conditions 
(A'^ = 10*, U — 0.69). Horizontal and vertical scales are as in 
Fig. I 

As the evolution proceeds, particles from the conden- 
sate "evaporate" into the gaseous phase 0, H. From 
Fig. H we note that, a,t t — 10'', only faint traces of the 
condensed phase remain in the ensemble. We have stud- 
ied the dynamics of this evaporation process as follows. 
First of all, we have numbered the rotators in each real- 
ization of the HMF model in such a way that, at the 
initial time, Pi{0) < Pj{0) for i < j. At subsequent 
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FIG. 7: Evolution of (a) the size of the condensate phase 
and (b) the average energy per particle in the condensate, 
for a single realization of the HMF model with A'^ = lO** and 
U = 0.69. Each dot corresponds to a running average over 
100 time units. The straight line in (a) has slope —0.85. 

The number n of rotators in the condensed phase dis- 
plays considerable oscillations in time. According to the 
discussion in the previous paragraph, we choose to mea- 
sure n in a single (typical) realization for TV = 10^ and to 
average over time intervals. Figure |^(a) shows the evo- 
lution of n averaged each 100 time units. The conden- 
sate size decreases steadily, displaying most of the time a 
power-law decay with a nontrivial exponent, n ^ t^'^'^^. 
This rather smooth decay is abruptly interrupted by a 
sudden collapse of the condensate, at t ^ 1.6 x IC*. In 
Fig. 0(b), we show the evolution of the average energy 
per particle in the condensate. We see that it varies non- 
monotonically with time, though its value is never far 
from u = 1. Other realizations for the same value of iV 
confirm that these results are generic. Moreover, they are 
robust with respect to the definition of the threshold S. 
For different values of N, realizations in the short-time 
range show that the collapse of the condensed phase oc- 
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curs at a time preceding, but approximately equal to, 
i,„i,i. In other words, when the minimum of the nonequi- 
librium temperature is attained, the condensate has just 
evaporated. From then on, the system enters the final 
stage of its relaxation to canonical equilibrium. 

In the long-time range, when the ensemble of rota- 
tors has presumably attained an ergodic state, its //-space 
configuration can be described statistically, in terms of 
averages over realizations, by means of a distribution 
function f{9,p)- Previous results on the //-space distri- 
bution of the HMF model refer, in all cases, to reduced 
distribution functions, feiO) and Jpip)- In Ref. ||, the 
momentum distribution function is used to illustrate, for 
U = 0.69, slow relaxation for long times. For other 
values of U, far from the critical point, both fo and 
fp were compared at moderately long times with the 
prediction of canonical statistics, showing the relatively 
fast attainment of canonical equilibrium for such ener- 
gies. In Refs. ^ and in contrast, deeper significance 
is attributed to nonequilibrium distribution functions — 
always in connection with QSS. There, LRT suggest that 
fp exhibits the fingerprints of nonextensive thermody- 
namics, and support this claim showing that the momen- 
tum distribution can be acceptably fitted, at low values 
of p, by the power-law distribution derived from Tsal- 
lis's nonextensive formulation Tsallis's formulation, 
however, predicts the onc-particlc probability distribu- 
tion for the energy, and not for the momentum. Now, 
since the one-particle energy |lf| 
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1 — Mx sin 6 — My cos ( 



(4) 



depends both on 9 and p, statistical descriptions in terms 
of e and p are equivalent only if the //-space distribution 
is independent of 9, i.e. if the energy distribution func- 
tion /e(e) can be unambiguously given in terms of fp{p)- 
This is not our case for U < Uc, unfortunately, since 
the angle distribution is not uniform, as illustrated by 
our Figs. 1^ and ^, and by LRT's figure 3(a) in Ref. |^. 
Ignoring the angular dependence of the one-particle dis- 
tribution in //-space invalidates the identification of fp 
as an equilibrium distribution derived from a generalized 
thermodynamical formulation like Tsallis's, which — just 
as ordinary thermodynamics — yields canonical distribu- 
tions as functions of the energy. 

We have evaluated by means of numerical realiza- 
tions of the HMF model for N — 10^. From the energy 
distribution function we can in turn calculate the energy 
probability density -P(e), since /^(e) = p{e)P{e), where 
p(e) is the density of states. The HMF density of states 
can be expressed in terms of elliptic functions, as 
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for 1 — M<e<l-|- M, where K and F are the com- 
plete and incomplete elliptic integrals of the first kind, 
respectively The energy is bounded from below, at 
Emin = 1 — M. Figure ^ shows the results for P(e) at three 
times, plotted as a function of e — emin- The straight line 
represents the canonical prediction, Poq(e) oc exp(— e/T). 
Observe that t = 5 x 10'* corresponds to the minimum 
of the nonequilibrium temperature, at the very begin- 
ning of the long-time range. For this time we note a 
considerable deviation from canonical equilibrium, with 
an abrupt cutoff at e — e„iin w 1 preceded by an over- 
populated interval. For very low energies, however, the 
coincidence with the equilibrium distribution is already 
quite good. For t = 2 x 10^ the cutoff persists but it is 
less abrupt and has shifted to higher energies, while in 
the low-energy range the distribution is closer to canon- 
ical equilibrium. Finally, for t — 10^, canonical equi- 
librium is well established in almost the whole range of 
energies shown in the figure. Deviations of the nonequi- 
librium temperature with respect to Teq are therefore to 
be ascribed to the depleted tail at high energies. From 
these results, we conclude that in the long-time regime 
the equilibrium energy distribution of the HMF model 
with WB initial conditions is built up from low ener- 
gies, with a depleted high-energy tail that recedes as time 
elapses. 
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FIG. 8: Probability distribution of energies at three times in 
the long-time range, for a system of A'^ = 10* rotators at [/ = 
0.69, as a function of e — tmin. The straight line stands for the 
exponential canonical prediction with temperature T = 0.476. 

We remark that LRT have assigned a distribution func- 
tion to states which are clearly in the short-time range. 
In that range the system is still strongly nonergodic and, 
as discussed above, such statistical description is there- 
fore of very limited significance. Their best fitting, for 
N — 10^ a.t t — 1200, corresponds in fact to a state 
where probably most of the ensemble is still in the con- 
densed phase. Moreover, they compare results for differ- 
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ent system sizes at the same time [see their figure 2(c) 
in Ref. ^, disregarding the scahng laws of time scales 
with the size, which threatens the validity of their evalu- 
ation of the fitting quality as a function of N [their figure 
2(d)]. Finally, we mention that the abrupt drop of fp for 
large momenta — also seen at larger times in our results 
for P(e), Fig. |— eludes any description based on Tsal- 
lis's thermodynamics. The artificial cutoff imposed by 
LRT to their fitting of /p, aimed at mimicking its abrupt 
drop, can hardly find any consistent justification within 
such formalism. We conclude that no evidence supports 
the association of noncquilibrium HMF states with gen- 
eralized thermodynamical equilibria. 

In summary, we have reexamined the slow-relaxation 
dynamics of the Hamiltonian mean-field model, when the 
evolution starts from water-bag initial conditions. First, 
we have revised the evidence that would support the exis- 
tence of quasistationary states in the noncquilibrium dy- 
namics of the model. In that regard, we have shown that 
the plateaus in the temporal evolution of the noncquilib- 
rium temperature T{t) [^0 are a mere artifact of the 
use of a logarithmic scale for the time axis in the plot of 
T{t). Moreover, a detailed analysis of the evolution of 
the nonequilibrium temperature shows that the dynam- 
ics can be decomposed into two well-defined dynamical 
regimes: a short-time range characterized by the pres- 
ence of a condensed phase, in which the system preserves 
strong memory of the initial conditions, and a long-time 
range, in which the condensed phase has evaporated and 
the system slowly relaxes toward equilibrium. A nat- 
ural boundary between these two regimes is given by 
the minimum reached by the nonequilibrium tempera- 
ture, which recedes toward larger times as the system size 
grows. The numerical evaluation of the probability den- 
sity function of the one-particle energy in the long-time 
range demonstrates that the system is slowly relaxing 
to the Boltzmann distribution, and does not support any 
connection with the generalized equilibrium distributions 
derived from Tsallis's nonextensive formalism. 



It must be stressed that neither previous work nor our 
present results elucidate the nature of the HMF evolution 
in the thermodynamical limit N —i- oo. The existence of 
different dynamical regimes, with different scaling laws 
in terms of N, makes an extrapolation from the available 
numerical results extremely uncertain. The possibility of 
performing significant, extensive numerical realizations 
of this model with current computational facilities seems 
to be restricted to ensembles of, at most, N ^ 10*^ el- 
ements. Numerical realizations for much larger systems 
may be necessary in order to attain a meaningful ap- 
proach to the thermodynamical limit. 

The origin of slow dynamics in the Hamiltonian mean- 
field model, also revealed by the study of aging phenom- 
ena and strong memory effects in the long-time range 
remains to be explained as well. A plausible argu- 
ment [^ may perhaps arise from the study of the role 
of the state of zero magnetization, which for U < Uc 
is thermodynamically unstable, and its connection with 
critical slowing down [^ — a connection which has often 
been invoked but not yet investigated. Revealing the 
basic mechanisms of slow relaxation in the Hamiltonian 
mean-field model should definitely clarify its relation to 
systems with some degree of structural disorder or frus- 
tration, such as glasses, that also exhibit slow-dynamics 
features. In turn, this would contribute to a unified view 
on nonequilibrium phenomena in systems with many de- 
grees of freedom. 



Acknowledgments 

We thank S. Ruffo for his motivating remarks on the 
HMF model, and C. Tsallis for comments on his own 
work. D. H. Z. acknowledges hospitality at the Abdus 
Salam International Centre for Theoretical Physics (Tri- 
este, Italy), where this work was partially carried out. 



[1] M. Antoni and S. Ruffo, Phys. Rev. E 52, 2361 (1995). [10] 
[2] V. Latora, A. Rapisarda and S. Ruffo, Phys. Rev. Lett. [11] 

80, 692 (1998). [12] 
[3] V. Latora, A. Rapisarda and S. Ruffo, Physica D 131, 

38 (1999). [13] 
[4] V. Latora, A. Rapisarda and S. Ruffo, Phys. Rev. Lett. 

83, 2104 (1999). 
[5] F. Tamarit and C. Anteneodo, Phys. Rev. Lett. 84, 208 

(2000). 

[6] V. Latora, A. Rapisarda and C. Tsallis, Phys. Rev. E 64, [14] 
056134 (2001). 

[7] V. Latora, A. Rapisarda and C. Tsallis, Physica A 305, [15] 
129 (2002). 

[8] M. Antoni, H. Hinrichsen and S. Ruffo, Chaos, Solitons 

and Fractals 13, 393 (2002). [16] 

[9] 



M. Montemu rro. F. Tamarit and C. Anteneodo, cond- 
mat/0205355i 



E 



C. Tsallis, Physica A 302, 187 (2001). 
C. Tsallis, Physica A 221, 277 (1995). 
T. Tsuchiya, T. Konishi, and N. Gouda, Phys. Rev 
50, 2607 (1994). 

Note that the definitions of the water-bag initial condi- 
tion respectively used for gravitational systems [|l^ and 
for the HMF model are substantially different. To our 
knowledge, the possible consequences of this difference 
have not been discussed in the literature. 
Lj. Milanovic, H. A. Posch, and W. Thirring, Phys. Rev. 
E 57, 2763 (1998). 

The analysis of S{t) — like most studies of the HMF 
dynamics — disregards the short initial stage of kinetic 
relaxation. 

Data obtained from the version of Ref. |6| at the public e- 



print archive http://arXiv.org/, cond-mat/010354C . The 
original figure also includes data for N — W, but they 



8 



do not span a sufficiently wide temporal range. 
For instance, the function T{t) = Ti + (To — 
Ti) exp(— t/r), whose derivatives are all monotonic func- 
tions, acquires an inflection point at f = r if plotted as a 
function of log t. 



[18] Data obtained from the version o f Ref. at the pub lic 
e-print archive http:/ /arXiv.org/ , |cond-niat/0109056 . 

[19] M. Abramowitz and I. Stegun, Handbook of Mathematical 
Functions (Dover, New York, 1972). 



